Saltation transport on Mars 
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We present the first calculation of saltation transport and dune formation on Mars and compare 
it to real dunes. We find that the rate at which grains are entrained into saltation on Mars is one 
order of magnitude higher than on Earth. With this fundamental novel ingredient, we reproduce 
the size and different shapes of Mars dunes, and give an estimate for the wind velocity on Mars. 
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The surprising discovery of sand dunes on Mars by 
Mariner 9 in 1971 has challenged geologists and plan- 
etary physicists over years vifith the following enigma: 
Could Mars dunes have been formed by today's martian 
thin atmosphere? Bagnold flj first noticed that the sand 
of aeolian dunes is transported through saltation, which 
consists of grains travelling in a sequence of ballistic tra- 
jectories and producing a splash of new ejected grains 
when colliding back onto the soil. But saltation occurs 
only if the wind friction speed — which, together with 
the air density pfluid, is used to define the air shear stress 
T = pfluidu^ — is larger than a minimal threshold u*t- 
Because the atmospheric density of Mars is almost 100 
times lower than on earth, only winds of strength one 
order of magnitude higher than that on our planet are 
capable of mobilizing the basaltic grains of mean diam- 
eter d = 500 /im 'i'l which constitute Mars dunes [3, 0] ■ 
Such winds occasionally can occur on Mars [1,[1,[3|- Here 
we find that the strong erosion due to these winds allows 
to understand dune formation on the red planet under 
present conditions. 

Once saltation starts, the number of saltating grains 
first increases exponentially due to the multiplicative pro- 
cess inherent in the splash events. But since the wind 
loses momentum to accelerate the grains, the flux satu- 
rates after a transient distance [1, 13, U , l2| • At sat- 
uration, the air shear stress Ta decreases to the threshold 
Tt ~ PfluidU*t, while the flux of grains increases to its 
maximum value. This transient phenomenon introduces 
the only relevant length scale in the physics of dunes: 
Any sand patch which is shorter than this saturation dis- 
tance will be eroded and disappear. On the other hand, 
large enough sand hills develop a slip face wherever the 
slope exceeds the angle of repose, 0^ = 34° [ij. 

Successful modeliiiK of the formation of dunes was re- 
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FIG. 1: Barchans in the Arkhangelsky crater, 41.0° S, 25.0° W 
on Mars: Mars Orbiter Camera (MOC) images on the left 
(image courtesy of NASA/JPL/MSSS) and calculated dunes 
on the right. Main plot: L vsW of the Arkhangelsky (circles) 
and north polar barchans at 77.6°N, 103. 6°W (stars). Calcu- 
lations of Arkhangelsky (north polar) dunes are represented 
by the continuous (dotted) line, obtained with ?i,/it*t ~ 1-45 
(1.80). The dashed line in the upper inset corresponds to ter- 
restrial dunes obtained with u,/u*t = 1.45. In the lower inset, 
we see L/Lmin vs H/Hmin from calculations of the Arkhangel- 
sky (triangles) and terrestrial dunes (dashed line). 



rate, T{x,y), which is the difference between the verti- 
cal flux of ejected grains and the vertical flux of grains 
impacting onto the bed: F = (f>{n — 1), where n is the av- 
erage number of splashed grains. At saturation (ra ~ Tt), 



can be written as 13|: 



cently achieved [IJ, The dune model consists of a 



system of continuum equations in two space dimensions 
which reproduce the shape of dunes, the wind profile and 
the sand flux and provide excellent quantitative agree- 
ment with measurements [l5l |. The sand bed can ex- 
change grains with the moving saltation layer. The local 
change in the flux, V-g, is proportional to the erosion 



n = 1 + 7( — 
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where 7, the entrainment rate of grains into saltation, 
determines how fast the system reaches saturation |13|. 
On the other hand, (f> = \q\/i, where £ is the average 
saltation length. The balance T{x,y) — V-(f yields a 
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differential equation for q: 
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where q^, = \q\{T — Tt) / (t — t^,) is the saturated flux [l3|, 
and i = vl^^{2vg/ g), where v°^° is the initial vertical ve- 
locity of the grains, Vg their mean velocity and g is grav- 
ity. Further, vf^ = aAvhon where At^hor is the horizontal 
velocity gain of the particle after one saltation trajectory, 
and a is an effective restitution coefficient for the grain- 
bed interaction In this manner, £ = {l/r)[2Vga/ g], 
with r = Vg/AvhoT, and the saturation length, 4, niay 
be then written as: 
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where 7 = r7. 

The dune model can be sketched as follows: (i) the 
wind shear stress over the topography is calculated using 
the algorithm of Weng et al. (1991) [1(5]. The wind ve- 
locity u over a flat ground increases logarithmically with 
the height z: u(z) = u^K~^\n.z/ zq, where k = 0.4 is the 
von Karman constant and zq is the aerodynamic rough- 
ness. A dune introduces a perturbation in the shear stress 
whose Fourier-transformed components are 
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where the coordinate axes x and y are parallel, respec- 
tively, perpendicular to the wind direction, and ky 



are wave numbers, = ^^k^ + k^ and e — 0.577216 

(Euler's constant). £ is the horizontal distance between 
the position of maximum height, i?max, and the posi- 
tion of the windward side where the height is iJniax/2 
[1^ . U{1) ^u{l)/u{hra) is the undisturbed wind veloc- 
ity at height / = 2K^C/\nl/zQ normalized by the veloc- 
ity at the reference height = C/ ^/logC/zg, which 
separates the middle and upper flow layers tl6]. The 
shear stress in the direction i {i = x, y) is then given by 
Ti^i [to(1 -I- Ti)], where tq is the undisturbed shear stress; 
(ii) next, the sand flux q (x,y) is calculated using eqs. 
(121) and ([3|); (iii) the change in surface height h{x,y) is 
computed from mass conservation: dh/dt — —V-q/psand, 
where Psand = 0.62pgrain is the bulk density of the bed 
(ist : and (iv) if sand deposition leads to slopes that lo- 
cally exceed the angle of repose, the unstable surface re- 
laxes through avalanches in the direction of the steepest 
descent, and the separation streamlines are introduced 
at the dune lee lj| . Each streamline is fitted by a third 



order polynomial connecting the brink with the ground 
at the reattachment point 1J|, and defining the "sepa- 
ration bubble", in which the wind and the flux are set 
to zero. The model is evaluated by performing steps i) 
through iv) computationally in a cyclic manner. 

In this letter, our aim is to reproduce the shape of Mars 
dunes using the martian present atmospheric conditions. 
Most of the relevant model parameters are known for 
Mars: The pressure, P, and temperature, T, which have 



been measured in several places on Mars 17| , determine 
Pfluid- Also known are the average grain diameter d and 
density, Pgrain — 3200 kg/m^, and gravity, g = 3.71 m/s^, 
while UH.t is calculated as in ref. [J]. 

One very common type of dune on Mars are harchans 
[1^ . They have one slip face and two horns, and propa- 

Ste on bedrock under conditions of uni-directional wind 
. Indeed, to develop a slip face, they must reach a 
minimum size. In the Arkhangelsky Crater on Mars (fig. 
1), the minimal dune width is Wmin ~ 200 m and the 
corresponding length Lmin ~ 400 m. On the other hand, 
Wrnin is around 13 times the saturation length, 4 [l9[. 
In this manner, Wmin yields, through eq. the wind 
strength u^,/u^,t on Mars. 

However, there is one unknown quantity for Mars 
which we need in order to solve the sand transport equa- 
tions: 7, which appears in eq. ([3]). Sauermann et al. 
(2001) [li] obtained 7 = 0.2 for sahation on Earth [H] by 
comparing with direct measurements of saturation tran- 
sients 0,1 111 . Such measurements are not available for 
Mars. 

We perform our calculations using open boundaries 
and a constant upwind u, in x direction at the inlet 
starting with a Gaussian hill, which evolves until dis- 
playing the linear relations between height H, width W 
and length L of barchans 13, 1^. 

For a constant 7, there is one value of which re- 
produces the minimal dune width Wmin — 200 m in the 
Arkhangelsky Crater, where P — 5.5 mb and T — 210 K 
17| . which gives u^t ~ 2.12 m/s. We obtained a surpris- 



ing result: If we take the same 7 as on Earth, the value 
of which gives Wmin — 200 m is around 6.0 m/s. This 
value is far too large to be realistic. The largestpeaks 
of martian range between 2.2 and 4.0 m/s 0,0]. To 
obtain the correct Wmin using m, within this maximum 
range, we must take 7 at least one order of magnitude 
larger than on Earth. We must find how much the mar- 
tian 7 = dn/d(ra/rt) (eq. Ilj) differs from the one on 
Earth. 

Anderson and Half (1988) [9] showed that the number 
of splashed grains is proportional to the velocity iiimp of 
the impacting grains. Let us rescale Uimp with Wcjo = 
\/'gd, which is the velocity necessary to escape from the 
sand bed Il2|. Further, Wimp scales with the mean grain 
velocity Vg [13[. In this manner, we obtain 7 oc Vgj \fgd. 

Moreover, Vg scales with [l3| and has only a very 
weak dependence on which we neglect, and thus 7 oc 
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Barchan field 


pfluid (kg/m^) 


u,t (m/s) 


(m/s) 


7 


Arkhangelsky 


0.014 


2.12 


17.8 


2.24 


77.6°N, 103.6°W 


0.022 


1.62 


12.3 


1.71 


Earth 


1.225 


0.22 


1.5 


0.20 



TABLE I: Main quantities controlling saltation on Mars 
and on Earth. it,t is around 80% of the direct entrain- 
ment threshold u,ft = A ^ Pgraingd/pfluid, where A — 

0. 129 [(1 + 6.0 X 10-Vpg.ain5d2-5)° V(l-928ReOtr - 1)°'] 

01, Re*ft = ittftcipfluid/'?, and rj is the viscosity of CO2 gas. 



u^t/Vgd- Since we know that 7 = 0.2 on Earth, where 
g = 9.81 ni/s^, d = 250 and M,t = 0.218 m/s [ll], we 
obtain 
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gd 
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Substituting eq. ^ into eq. ([3]), we obtain a closed set 
of sand transport equations from which the value of the 
wind friction speed on Mars can be determined. 

Eq. ^ gives 7 « 2.24 in the Arkhangelsky Crater, 
and Wmin = 200 m is reproduced for u*/u*t = 1-45 or 
u* = 3.07 m/s, which gives is ~ 15.5 m. Fig. 1 shows 
that — 3.07 m/s reproduces not only the minimal 
dune but also the dependence of the dune shape on its 
size. Furthermore, we tested the scaling relation ([6|) with 
a second martian barchan field which is near the north 
pole, and where Wmin ~ 80 m. At the location of the 
field, P = 8.0 mb and T = 190 K [13], and thus u*t « 
1.62 m/s. From Wmin, we obtain £s ~ 6.0 m, which gives 
w*/w*t ~ 1-8 or = 2.92 m/s. The plot in fig. 1 shows 
that the behaviour L against W of the barchans in both 
studied fields is well captured by the model. 

It is interesting that the obtained for the north polar 
field is very similar to that in the Arkhangelsky Crater, 
although is lower in the north polar field due to the 
higher pfluid (table 1). 

Summarizing, we discovered that the rate 7 at which 
grains enter saltation on Mars is about one order of mag- 
nitude higher than on earth. Taking this finding into ac- 
count, we found that the wind velocity which reproduces 
the size and shape of martian barchans is ~ 3.0 ± 0.1 
m/s, which is well within the range of estimated wind 
speeds on Mars d, Q . We found that Mars grains travel 
with a velocity Vg ten times higher than that of terres- 
trial grains, as observed in wind tunnel experiments for 
Mars [20|. Such high-speed grains produce much larger 
splash events than on Earth 2l|, and lead to higher 7 
values (table 1). 

Our results show that the larger splash events have a 
crucial implication for the formation of dunes on Mars. 
While on one hand the lower martian g and pfluid result 
in longer grain trajectories [20|, the higher rate at which 
grains enter saltation on Mars shortens the saturation 
transient of the flux, which determines the minimal dune 
size. The distance of flux saturation, As, is around six 
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FIG. 2: Dune velocity v as function of dune length L. We 
see that Mars dunes (filled symbols) move typically ten times 
faster than Earth dunes (empty symbols) of same L, obtained 
with similar values of u*/u*t as on Mars. 



times the characteristic length is [13|, which effectively 
increases with i. Here we found that As on Mars is indeed 
reduced by the larger martian splash. This explains the 
previously reported [22] failure of the scaling Wmin oc i 
for Mars. 

The values of M*/u*t obtained for Mars are within the 
range of the ones measured in terrestrial barchan fields 
[lit . Indeed, we see in fig. [5] that, for the same value of 
u*/u*tj Mars barchans would move ten times faster than 
those on Earth. 

Many exotic and up to now unexplained dune forms 
have also been observed on Mars. Dunes as those in 
fig. 3 cannot be obtained from calculations with a uni- 
directional wind. We found that the dune shapes in fig. 
3 may be obtained with a bimodal wind regime. In our 
calculations, the wind alternates its direction periodically 
with frequency 1/Tw forming an angle 2\9^\ as sketched 
in fig. 3a'. In both directions the strength is the same, 
namely = 3.0 m/s. In this manner, the value of u*/u*t 
is particular to each field, since M*t depends on the field 
location (table 2). To simulate the change of wind di- 
rection, we rotate the field by an angle 29^, keeping the 
wind direction constant. The separation bubble, thus, 
adapts to the wind direction after rotation of the field. 
We use open boundaries as in the calculations of barchan 
dunes. Initial condition is a gaussian hill as before, whose 
volume is taken according to the volume of the dune. 

The angle 20^ between the wind directions determines 
which of the different forms in fig. 3 is obtained. We 
found that a barchan moving in the resulting wind direc- 
tion is always obtained if ]6'w] < 50°. li \9^\ 50°, the 
dune shape in fig. 3a' is achieved. And for \9^\ between 
50° and 70°, elongated dune forms as those in fig. 3b' 
are obtained, which elongate further and further in time. 
As \9v,\ — > 90°, a barchanoidal form of alternating slip 
face position appears. 

Moreover, we found that the structure observed in the 
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FIG. 3: MOC images on top (courtesy of NASA/JPL/MSSS) 
and calcuiations obtained witli bimodal wind regimes on bot- 
tom. A sketcfi sliowing tfie definition of tfie angie 6^ of tlie 
wind direction (arrows) is sliown in a'. Tlie wind clianges 
its direction witli frequency l/Tw. We cliose Tw = 2.9 days, 
5.8 days and 0.7 day to obtain dunes in a', b' and c', respec- 
tively. Tfie dune in a' lias been obtained with |6wl = 50° and 
the dune in b' with |Sw| = 70°. In c', a linear dune obtained 
with |0w| = 60° decays into barchans after the angle \6^\ is 
reduced to 40° . Dune in b' elongates to the right further and 
further with time, which is not observed for dune in a'. Dune 
in c' decays further until only a string of rounded barchans is 
visible. 



Field 


location 


Pfluid (kg/m-*) 


Utt (m/s) 


fig. 3a 


48.6°S, 25.5°W 


0.017 


1.89 


fig. 3b 


49.6°S, 352.9°W 


0.014 


2.06 


fig. 3c 


76.4°N, 272.9°W 


0.03 


1.35 



TABLE II: For each dune field in fig. 3, the fluid density pfluid 
and the threshold u,t is calculated from the local pressure and 
temperature which are taken from ref. In spite of the 

broad range of u,t, all dune forms in flg. 3 have been obtained 
with one single value of u« = 3.0 m/s. 

dune field of fig. 3c can be obtained by a change in the 
local wind regime. The dune shape in fig. 3c' has been 
obtained in the following manner: (i) first, an elongated 
dune form as the one in fig. 3b' is formed with an angle 
\6yj\ = 60°; (ii) next, the angle |6'w| has been suddenly 
reduced to 40°. Thereafter, the linear dune becomes un- 
stable and decays into a string of rounded barchans seen 
in fig. 3c. 

Each one of the dune shapes in fig. 3 can only be 
achieved if the time Tw is of the order of 1 — 5 days. If 
the period is too large, of the order of 1 month, then the 
dune evolves into a barchanoidal form. 

In conclusion, we found from calculations of martian 
barchans using the present atmospheric conditions of 
Mars, that the splash on Mars must be about ten times 



larger than on Earth. Furthermore, we obtained a gen- 
eral equation for the entrainment rate of grains into salta- 
tion which could be used in future calculations of other 
dune forms on Mars, Venus or Titan. We also found that 
winds on Mars don't exceed 3.0 ±0.1 m/s, which can ex- 
plain the shape of elongated dunes formed by a bimodal 
wind. It would be interesting to make a full microscopic 
simulation for the saltation mechanism of Mars similar 
to the one that was recently performed by Almeida et al. 
(2006) [l^l to confirm our findings microscopically. 
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